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ABSTRACT 

There have now been three supernova-associated 7-ray bursts (GRBs) at redshift 
z < 0.17, namely 980425, 030329, and 031203, but the nearby and under-luminous 
GRBs 980425 and 031203 are distinctly different from the 'classical' or standard GRBs. 
It has been suggested that they could be classical GRBs observed away from their jet 
axes, or they might belong to a population of under-energetic GRBs. Recent radio 
observations of the afterglow of GRB 980425 suggest that different engines may be 
responsible for the observed diversity of cosmic explosions. Given this assumption, 
a crude constraint on a luminosity function for faint GRBs with a mean luminosity 
similar to that of GRB 980425 and an upper limit on the rate density of 980425- 
type events, we simulate the redshift distribution of under-luminous GRBs assuming 
BATSE and Swift sensitivities. A local rate density of about 0.6% of the local super- 
nova Type Ib/c rate yields simulated probabilities for under-luminous events to occur 
at rates comparable to the BATSE GRB low-redshift distribution. In this scenario the 
probability of BATSE/HETE detecting at least one GRB at z < 0.05 is 0.78 over 4.5 
years, a result that is comparable with observation. Swift has the potential to detect 
1-5 under-luminous GRBs during one year of observation. 
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Discovery of GRB 980425, associated with the very nearby 
supernova (SN) 1998bw (at z — 0.0085, correspond- 
ing to about 40 Mpc), heralde d a new era in und er- 
standing the origin of GRBs iGalama et alJ Il998allbh . 
Recently, the detection of spectroscopic features in the 
light curve of GRB 03032 9 , similar to those see n 
in SN 1998bw iStanek et alJ 120031: iHiorth et all 12003!) . 
has strengthened the SN 1998bw/GRB 980425 associa- 
tion. These observations support th e 'collapsar' model 
jZhang. Wooslev fc MacFadvenl 120031) in which a Wolf- 
Rayet progenitor, possibly in a binary system, undergoes 
core collapse, producing a compact object surrounded by 
an accretion disk, which injects energy into the system and 
thus acts as a 'central engine'. The energy extracted from 
the system gives rise to a Type Ib/c SN explosion and drives 
collimated jets along the progenitor rotation axis, produc- 
ing a prompt GRB and afterglow emission-see the review 
by Zhang & Meszaros (2004). 

Discovery of the GRB-SN association was an impor- 
tant breakthrough, but GRB 980425 had an unusually low 
luminosity — it was under-luminous in 7-rays by three orders 
of magnitude compared to 'classical' GRBs. It was suggested 
that it could be a very rare event and not a member of the 
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classical GRB population. This explanation seems unlikely 
given the discovery of the under-luminous and nearby GRB 
031203 (z = 0.105), the first analogue of GRB 980425. 

There are presently three GRBs (980425, 030329, and 
031203) de finitely associated with extremely energe tic Type 
Ib/c SNe JProchaskal 120041: iMalesani etatl I 20041) . all oc- 
curring at z < 0.17. GRB 030329 is classified as classi- 
cal in the context of energy emission, but was relatively 
close at z = 0.17. However, it is difficult to reconcile the 
under-luminous GRBs 980425 and 031203 with the classical 
population. One simple explanation is that under-luminous 
bursts are GRBs observed away from the jet axis. 



iGuetta et all J2004I) argue that a unified picture can 
only be obtained by using a luminosity function (LF) that in- 
cludes all luminosities down to that of GRB 980425, so that 
the probability of observing the three low-z events is non- 
negligible. They show that for GRBs 980425, 030329, and 
031203 to belong to the classical burst population, the LF 
must be a broken power-law. This is an attractive proposal 
in that GRBs 980425 and 031203 can be explained within 
the bounds of currently popular GRB progenitor models by 
extending the LF to accommodate GRB 980425. They cal- 
culate that if this is the case, no bright burst within z — 0.17 
should be observed by a HETE-like instrument within the 
next ~ 20 yr. 
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2 ANOMALOUS GRBS 

2.1 Evidence for intrinsically sub-energetic events 

ISoderberg et, alJ lj2004) argue that if GRB 031203 was ob- 
served 'off axis', then the radio afterglow should brighten 
as the ejecta slows down, but they did not observe any re- 
brightening. They find that the afterglow is faint, indicat- 
ing that the explosion was under-energetic. Similarly, there 
is no evidence of re-brighteni ng for GRB 980 4 25, d espite 
radio calorimetry since 1998. ISoderberg et alJ fc0Q4ft con- 
clude that GRBs 980425 and 031203 were intrinsically sub- 
energetic events. 

If SN 1998bw was a rare and un usually sub-energ e tic SN 
distinct from local SNe and GRBs, ISoderberg et alJ (120041) 
claim the characteristics of SN 1998bw/GRB 980425 are not 
a result of the observer's viewing angle but of the properties 
of its ce ntral engine. SN 199 8bw was an engine-driven ex- 
plosion (iLi fc Chevalierl lT999). in which 99.5% of the kinetic 
energy (~ 10 50 e rgs) was coupled to relativistic ejecta of 
Lorentz factor 2 llKulkarni et al] IT998). while a mere 0.5% 
went into the ultra-relativistic flow. In contr ast, 'classical' 
GRBs c ouple most of their energy into 7-rays. iBerger et alJ 
J2003bl) claim that the observed diversity of cosmic explo- 
sions (SNe, x-ray flashes (XRF), and GRBs) could be ex- 
plained with a standard energy source but with a varying 
fraction of that energy injected into relativistic ejecta. Dif- 
ferent engines may be responsible for the observed diversity 
of cosmic explosions, implying that classical GRBs represent 
one class of event, one in which 7-rays channel most of the 
energy away from a central engine. 

It is evident that SN 1998bw could be a member of a 
distinct class of SN explosions. But how rare is SN 1998bw 
in the context of T ype Ib/c SNe and classical GRBs? 
IBerger et all j2003al) carried out a systematic program of 
radio observations of Type Ib/c SNe using the Very Large 
Array to place the first constraint on the rate density of SN 
1998bw type events. Of the 33 SNe observed from late 1999 
to the end of 2002, they conclude that the fraction of events 
similar to SN 1998bw is at most 3%. Furthermore they find, 
by comparison of the SN radio emission to that of GRB af- 
terglows, that none of the observed SNe could have resulted 
from a classical GRB. 

2.2 Evidence for an off-axis model for GRB 
031203 

The evidence that GRB 031203 was an intrinsically faint 
and nearly spherical exp losion is not widely accepted. 
iRamirez-Ruiz et al.l (120041) disagree with this interpretation 
and argue, from models fitted to the observed x-ray light 
curve and the radio afterglow, that GRB 031203 was a classi- 
cal GRB viewed off axis. They find that most spherical mod- 
els under-predict the x-ray flux at late times by at least two 
orders of magnitude, and prefer to interpret GRB 031203 as 
a highly collimated GRB viewed off-axis. 

For the case of an observer located outside the jet aper- 
ture, # bs > do, the prompt GRB emission and its early af- 
terglow are considerably weaker than for on axis, # b s < do- 
An observer at dobs > do sees a rising afterglow light curve 
at early times, which approaches that seen by an on-axis 
observer at late times. The emission remains low until the 
cone of the beam intersects the observer's line of sight. 



In the off-axis jet scenario, with viewing angle 6* bs ~ 
2#o, GRB 031203 can be model led as a GRB viewed a few 
degrees outside of a conical jet. IRamirez-Ruiz et alJ (I2004T) 
show that if GRB 031203 is modelled as an off-axis obser- 
vation, its energy emission in 7-rays is about 10 53 ergs, con- 
sistent with classical GRBs. 

The off-axis model of IRamirez-Ruiz et alJ i2004l) used 
to explain the weak afterglow of GRB 031203 casts some 
doubt on the claim that it is an intrinsically weak event. Fur- 
thermore, because GRBs 031203 and 980425 share a com- 
mon deficit in 7-ray emission, it is possible that these 'out- 
liers' are a result of off-axis observations. The possibility 
that GRBs may represent a class of cosmic explosions with 
a broad range of energies provides the basis for the following 
simulation. 



2.3 A population of under-luminous GRBs 

Given the potential of the Swift satellite, a multi-wavelength 
GRB observatory (http://swift.gsfc.nasa.gov/) launched on 
2004 November 20, to localize hundreds of GRBs, observa- 
tions of events similar to GRB 980425 will provide new in- 
sight into GRB progenitor populations. This is strong moti- 
vation to constrain the probability of detecting similar GRB- 
SNe assuming a simple, if highly uncertain, model. The first 
upper limits on the rate density of Type Ic SNe associated 
with GRBs are being made, based on the presence of jetted 
emissions (e.g. Berger et al. 2003a). Given such upper limits, 
one can at least provide further constraints on GRB pro- 
genitor populations using continued satellite observations. 
Furthermore, we provide an example to demonstrate how 
an under-luminous GRB population would manifest during 
the BATSE observation period and the era of Swift. 

Firstly, we assume GRB 980425 is a 'typical' member of 
a class of relatively rare GRB (compared to classical GRBs) . 
Based on the luminosity of GRB 980425, the mean luminos- 
ity of this population is about 3 orders of magnitude less 
than that of the classical population. Secondly, the local 
rate density must be less than 3% of the Type Ib/c SN rate. 
Another observational constraint is based on the observed 
rate of GRB 980425 type events out to z — 0.0085 assuming 
a 4.5 yr (BATSE) observation period; the probability of oc- 
currence of this very nearby GRB must be compatible with 
the BATSE GRB distribution. Finally, we assume the 7- 
ray emissions are isotropic (not beamed) based on the radio 
observations of the aft erglows of GRBs 980425 and 031203 
iSoderberg et af .120041) . With these constraints, we simulate 
the observed GRB distribution for BATSE and Swift sensi- 
tivities. 



3 THE GRB LUMINOSITY FUNCTION 

The GRB LF, together with the flux sensitivity threshold 
of an instrument, determine the fraction of all GRBs poten- 
tially detectable with that instrument: 

r 00 

^grb(^) =S d / p(L)dL, (1) 

where V'grb(-z) is the GRB rate scaling function, Sd is the 
fraction of sky that the detector scans, and p(L) is the 
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GRB LF with L the intrinsic luminosity in units of pho- 
tons s~ . With fnm denoting the instrumental flux sensitiv- 
ity threshold, in photons s _1 m -2 , the minimum detectable 
luminosity can be expressed as a function of redshift by 
Lih-n(z) = 47r_DL 2 (z)/iim, with -Dl(z) the luminosity dis- 
tance. 

Most models for the clas sical GRB LF are ba sed on the 
luminosity-redshift relation iSchaefer et al-lfeoOll) . But the 
models are biased by the redshift sample and the sensitivity 
limit related to the redshift estimate. This limit has to be 
lowered at least by an order of magnitude to encompa ss the 
complete range of luminosities. iFirmani et alJ (120041) show 
that by jointly fitting to the observed differential peak-flux 
and redshift distributions, the best fit for the LF takes a 
form that evolves weakly with redshift. However, there is no 
consensus on the form of LF for classical GRBs; possibilities 
include a single power law, a double power law and a log- 
normal distribution. 

For an under-luminous population of GRBs modelled on 
the single GRB 980425, the choice of LF is so uncertain that 
the form of the function is somewhat arbitrary. Nonetheless, 
for definiteness and for comparison, we use log-normal distri- 
butions for the classical GRBs and an under-luminous popu- 
lation, both with observation-based statistical moments that 
fit the data: 

2 /2 



p(L) = 



V2I 



: exp 



ML/U)] 2 
2a 2 



J_ 
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where a and Lq are the width and average luminosity, re- 
spectively. We take a — 2 and L = 2 x f0 5B s _1 for the 
classical GRBs, with /u m = 0.2 and 0.04 photons s _1 cm~ 2 
for BATSE (classical G RBs) and Swif t respectively, and take 
about O.f for S d l|Guetta et alJl2004t l. Assuming that GRB 
980425 is representative of an under-luminous population, 
we take Lo = 2 x 10 53 s _1 , three orders of magnitude less 
than for the classical GRBs. 



4 GRB RATES AND DETECTION 
PROBABILITY 

One can express the differential GRB rate in the redshift 
shell 2 to z + dz as 



dR = 4>{z)— „ dz , 



(3) 



dz 1 + z 

where dV is the cosmology-dependent co-moving volume el- 
ement and R(z) is the GRB event rate, as observed in our 
local frame, for sources out to redshift z. Source rate density 
evolution is accounted for by the dimensionless evolution fac- 
tor e(z), which is normalized to unity in the present-epoch 
universe (z = 0), and ro is the z — rate density. The (1 + z) 
factor accounts for the time dilation of the observed rate by 
cosmic expansion. 

We assume a 'fiat-A' cosmology with Q m = 0.3 and 
SIa = 0.7 for the present-epoch density parameters, and 
take Ho = 70 km s _1 Mpc -1 for the Hubble parameter at 
z = 0. The star form ation rate (SFR) model SF2 of 
iPorciani fc Madaul ll200 if) is re-scaled to this cosmology and 
converted to a normalized evolution factor e(z). This SFR 
model levels off to an essentially constant rate density, of 
order 10 times the z — value, at z > 2; a star formation 
cutoff at z = 10 is assumed. 



For the classical GRBs we use ro = 0.9 yr" 1 Gpc" 3 , 
obtained by scaling the all-sky Universal rate to 692 yr -1 , 
the value implied by the GUSBAD 1 catalogue . This is com- 
parab le to the value of 1.1 yr -1 Gpc~ 3 from iGuetta et ail 
(2004), obtained using a different SFR and a broken power- 
law luminosity function. 

As GRBs are independent of each other, their distri- 
bution is a Poisson process in time: the probability for at 
least one event to occur in the volume out to redshift z dur- 
ing observation time T at a mean rate R(z) is given by an 
exponential distribution: 

-R(z)T 



p{n>l;R{z),T) = 1 



(4) 



This formula can be used to define a 'probability event 
horizon' — as observation time increases, how often will rarer, 
more local events, be observed? See ICoward fc Burmanl 
( 2005) for a description. Based on equation @, the probabil- 
ity of at least one GRB occurring in z < 0.17 during 4.5 yr is 
about 0.5, implying that the observed GRBs in this volume 
need not be considered anomalous. But the probability of a 
GRB occurring in z < 0.01 over 4.5 yr is 0.00015, implying 
that GRB 980425 (z = 0.0085) is either an extreme outlier 
or a member of a different GRB population. We model the 
GRB redshift distribution under the latter assumption, with 
the constraint that the probabilities and rates are consistent 
with the observed BATSE GRB distribution. 



5 SIMULATING AN UNDER-LUMINOUS GRB 
REDSHIFT DISTRIBUTION 

A GRB distribution comprised of two populations, with dif- 
ferent local rate densities and mean luminosities can be ex- 
pressed as: 



dR = 



dV e(z) 



dz 1 + z 



^c{z)r c +ip n (z)ro 



dz , 



(5) 



where tj) c (z)ro an d V , u( z )'"o are the scaling functions and lo- 
cal rate densities of the classical and under-luminous GRBs 
respectively, ft is assumed that both rates follow the SFR 
density so that e(z) is the same for both populations. 

Figure 1 plots ip c (z) and ip u (z) for BATSE and Swift 
sensitivities. It is evident that at z > 0.1, the flux-limit 
of the detectors severely limits the potential detectability 
of the under-luminous GRB population (mean luminosity 
Lo = 2 x 10 53 s _1 ). For z < 0.1 the scaling function is 
non-negligible even though events will be rare inside such 
a relatively small volume. The sensitivity of Swift implies 
that it could potentially detect most under- luminous events 
occurring inside a volume bounded by z = 0.01. 

Clearly the observed redshift distribution of GRBs and 
expected rates for very energetic Type Ib/c SNe do im- 
pose constraints on the local rate density rg. Importantly, 
GRBs 980425 and 031203 show no evidence for jets, im- 
plying that there is no geometric rate enhancement factor 
required to account for unseen bursts of similar type. As 
a first approximation to rj}, we take the classical rate r§ 
increased by the beaming rate enhance ment factor for clas - 
sical bursts, using a value of 250 from lFrail et al.| 1120011) . 
This gives r^ « 220 yr -1 Gpc -3 , which is about 0.6% 

1 http:/ /www. astro. caltech.edu/~mxs/grb/GUSBAD/ 
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Redshift, z 

Figure 1. GRB scaling functions ip c for the classical population, 
and ipu for an under-luminous population, shown using BATSE 
and Swift sensitivities. ip u is calculated using the assumption that 
GRB 980425 is a typical event from an under-luminous population 
of mean luminosity Lo = 2 X 10 53 s — 1 , three orders of magnitude 
less than the mean luminosity of the classical GRBs. A value of 
0.1 for the scaling function corresponds to all GRBs in the field 
of view of the detector (0.1 str) being potentially detectable. 



Table 1. The probability of observing at least one GRB inside 
volumes bounded by z = 0.0085, 0.05, and 0.17 during 4.5 yr 
of observation using a single GRB distribution (classical) with 
BATSE sensitivity and a double distribution (classical + under- 
luminous) with BATSE and Swift sensitivities, (labelled as Dou- 
ble (BATSE) and Double (Swift) respectively. 

p(z < 0.0085) p(z < 0.05) p(z < 0.17) 

Classical (BATSE) 0.00015 0.02 0.5 

Double (BATSE) 0.04 0.78 0.99 

Double (Swift) 0.04 0.96 0.99 



of the l ocal SN Type Ib/c rate, r^ N ' bc sa 3.7 x 10 4 yr" 1 
Gpc" 3 llzzard. Ramirez-Ruiz fc Toutll2004D — a result that 
supports the view that SN1998bw was a relatively rare 
and unusually energetic SN. We note that the quoted SN 
Type Ib/c rates may be underestimated because many core- 
collapse SN are lost to extinction in most surveys to date. 

The mean of the luminosity function needs to be re- 
duced by 3-4 orders of magnitude and the variance reduced 
from 2 to 1.5 to crudely fit the resulting probability distribu- 
tion with the observed rates of GRBs at small redshift; that 
is, the probability of occurrence of the very nearby GRB 
980425 (z = 0.0085) must be compatible with present ob- 
servations. For this condition to be satisfied we find that 
a must be smaller than that used for modelling the clas- 
sical population (a = 2), otherwise the rate of events at 
small z would be too large. If a broader luminosity distri- 
bution (larger variance) is employed, for example a = 2.5, 
the mean luminosity must be reduced to 4 orders of magni- 
tude less than the mean classical GRB luminosity, to yield 
observationally consistent probabilities. For definiteness and 
consistency, we assume a = 1.5 and Lo = 2 x 10 53 s _1 in all 
calculations 




Redshift, z 



Figure 2. The differential number of GRBs as a function of red- 
shift for an observation time of 1 yr using the three models de- 
scribed in Table 1. The BATSE classical and double distribution 
models predict similar numbers from z = 0.2 — 10. 




Figure 3. As for Figure 2 but plotting the cumulative number 
of GRBs as a function of redshift for an observation time of 1 
yr. A comparison of the double and classical models for BATSE 
shows that the under-luminous population causes a significant 
increase in numbers at z < 0.1, resulting in a higher probability 
of observing small-2 GRBs. Swift could potentially detect over 
100 GRBs during one year with up to 5 in z < 0.1 



Table 1, using the same parameters as figure 1, shows 
the probabilities for detecting at least one GRB in volumes 
bounded by z — 0.0085, 0.05, and 0.17 during 4.5 yrs ob- 
servation, assuming the sensitivities of both BATSE/HETE 
and Swift. Equation 10 is used to calculate the rates based 
on a classical distribution and equation © for a distribu- 
tion comprised of both classical and under-luminous bursts. 
The double distribution model increases the detection prob- 
ability of the very nearby GRB 980425 (z = 0.0085) to a 
still small (0.04), but significant, level compared with the 
extremely small probability (0.00015) from the classical dis- 
tribution alone. 
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Figures 2 and 3 plot the number distribution and cu- 
mulative number of GRBs observed over a 1 yr period for 
the classical and double distributions assuming BATSE and 
Swift sensitivities. For BATSE sensitivities it is evident that 
the under-luminous GRBs have no effect on the observed 
distribution at z > 0.2. They do contribute to the cumu- 
lative number at z < 0.2, a result that is compatible with 
the detection of GRBs 980425 and 031203. The cumulative 
number increases from about 1 at z = 0.1, for a BATSE 
sensitivity, to about 4 or 5 for a Swift sensitivity. 



6 CONCLUSIONS 

We have shown that a population comprised of classical and 
under-luminous GRBs is compatible with the currently ob- 
served GRB redshift distribution that includes the nearby 
and under-luminous GRBs 980425 and 031203. We find that 
a local rate density for an under-luminous GRB population 
of rg ~ 220 yr -1 Gpc" 3 , which is about 0.6% of the local SN 
Type Ib/c rate, rQ NIb/c w 6.3 x 10 4 yr" 1 Gpc" 3 , fits the ob- 
served low-redshift GRB distribution. Assuming that GRB 
980425 is typical in luminosity for an under-luminous popu- 
lation, we make a first crude constraint on such a population 
by taking a mean luminosity of Lo = 2 x 10 53 s _1 — about 3- 
4 orders of magnitude less than the mean luminosity of the 
classical GRBs. These two constraints yield a probability of 
BATSE/HETE detecting at least one GRB at z < 0.05 to 
be 0.78 over 4.5 years, a result that is compatible with the 
presently observed low-redshift GRB distribution. 

GRBs 980425 and 031203 may only appear faint 
and anomalous because of the sensitivity limit of 
BATSE/HETE. If such an under-luminous population is 
present, the increased sensitivity of Swift should enable it 
to detect and localize 5 times more GRBs than BATSE at 
redshifts z < 0.1. 

It seems reasonable that the observed broad distribu- 
tion of observed GRB luminosities may represent related, 
but different classes of engine-driven emissions powered by 
rotating massive compact stellar remnants. The definite as- 
sociation of some nearby GRBs with Type Ib/c SNe-a SN 
type that exhibits considerable diversity-supports the idea 
that there could be a diverse class of inner engines driving 
at least a fraction of GRBs. These classes may even form 
a continuum that encompasses Type Ib/c SNe, XRF, faint 
GRBs and classical GRBs. Hence the distribution of sources 
in redshift in this scenario would consist of the sum of the 
individual populations with different inner engines and lo- 
cal rate densities. There is most likely overlap between the 
emission characteristics of the various sub-populations, so 
much so that they may form a continuum of luminosities 
ranging from XRF to hard GRBs. 

If the anomalous GRBs are all shown to be off-axis ob- 
servations of classical GRBs, then the evidence for different 
classes of cosmic explosions related to GRBs will become 
more tenuous. Swift should provide new information on the 
diversity of inner engines and on the progenitors that pro- 
duce them, providing a wealth of data to help solve the 
cosmic riddle of identifying GRB progenitor populations. 
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